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Introduction 



* High speed aircraft are highly desirable for military and 
commercial applications. 

• “Hypersonic” speed is frequently defined as Mach No. > 5 


SR71 (Mach 3+): 


XI 5 (Mach 6.7): 




Introduction 





• Hypersonic air-breathing propulsion is one of the most difficult 
R&D challenges facing the worldwide aeronautics community. 

X51 Waver ider: 



• 1 out of 3 flights resulted in sustained combustion. 

• Hydrocarbon fueled. 
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Introduction 




• Current R&D uses ground tests and computational fluid 
dynamics (CFD) to prepare for flight tests. 


X43: 



• Flight 2 achieved Mach 7; Flight 3 achieved ~ Mach 10. 

• Hydrogen fueled. 
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Introduction 
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• High Speed / Hypersonic R&D uses flight tests, ground tests, 
and CFD - — all have difficulties: 

1. Flight tests: Expensive, difficult to instrument, very harsh 
environments for experimental aircraft that naturally lead to high 
failure rates. 

2. Ground tests: Difficult to achieve actual flight conditions, facility 
interference effects, difficult to instrument. 

3. CFD: Physical modeling is highly unproven - at best. Modeling 
of turbulence, chemistry, heat transfer, fluid-thermal-structural 
interactions is subject to high uncertainty. 
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Introduction 
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Turbulence models remains one of the key pacing technologies 
in Computational Fluid Dynamics (CFD). 


• An overview of key turbulence modeling areas for propulsion 
flows is presented. 

• Emphasis is placed on “practical” state-of-the-art today: 

- Standard practices using primarily RANS. 

- Promising new technology (i.e. LES, hybrid RANS/LES) that may be available for 
production use in near future. 

- Key shortfalls for which R&D is necessary. 

• Focus is placed on high-speed propulsion systems (i.e. 
scramjets); turbine engines are also addressed in less detail. 
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Key Turbulent Features of 
Scramjet Flowpaths 



HSST 




FOREBODY/INLET: 

Laminar-to-turbulent 
transition, Shock wave / 
turbulent boundary layer 
interactions 


ISOLATOR: Shock wave/ 
turbulent boundary layer 
interactions 


COMBUSTOR AND 
EXPANSION SYSTEM: 

3D compressible mixing, 
reacting flow, turbulent / 
chemistry interactions 
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Key Turbulent Features of 
Turbine Engine Flowpaths 



INLET: 

Transition, 

Separation 



COMPRESSOR: 
Swirling 3D flow, 
wakes, shock- 
interactions 


TURBINE: 
COMBUSTOR: Transition, 3D, very 

3D reacting flow, hi 9 h heat transfer, 
turbulent / chemistry tilm cooling 
interactions, multi-phase 


NOZZLE/MIXER, 

PLUME 

3D Turbulent Mixing, 
Compressibility, 
Acoustics 


8 



Presentation Outline 
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• Overview of Turbulence Modeling in Use for Propulsion Flows 

- RANS 

- LESand DNS 

• Boundary Layer Transition - Inlets and Turbines 

• 3D Boundary Layer Effects 

• Shock-Wave /Turbulent Boundary Layer Interactions 

• Combustor / Reacting Flows 

- Scalar Transport 

- Turbulent / Chemistry Interactions 

• Exhaust System Modeling 

- Jet and Mixing - RANS 

- LES-based Methods 

• Experimental Validation Data Needs 

• Conclusions 
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Navier-Stokes Equations 




dp 

dt 


+ 


oUi) = 0 
03C i 



dP dz 
— + — L 
dx i dXj 


Continuity 


Momentum 


dE d d dqj 

-^ + —-( u i ( E , + p )) = ^(u,T iJ )-- 1 - Energy 

dt d)Cj dXj dxcj 


r 1 
E, =pe + -pu i u i 
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Navier-Stokes for Perfect Gas 




P = pRT 


Ideal Gas 


2 dUj 


Newtonian Fluid 


s -I 

" 2 


3a, du j 
— L + — - 

\&j j 


Rate of Strain 


q j =-k 


dT_ 


Fourier’ s Law 
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Reynolds Averagi 




Traditional 

Reynolds-averaging 

Density weighting 
(Favre averaging) 


(Favre averaging) 
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Modeled Turbulent Terms - RANS 
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v 

-pii;u'. - t 
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Reynolds Stress 


-pu',h , = q T , 


Turbulent heat flux 


These terms replace ALL turbulent effects in the momentum and energy 

equations. 


Spatial Filtering 


E\ 1 1 ( 





^AOTV.C^ V 


/(-1,0 = / D C(* 


A)f(%,t)d 


Filtering operation 



G(x-^A)d i % = \ 


Filter function 



(Favre filtering) 


Note that the Favre spatial filtering here and Favre Reynolds averaging 
look similar but refer to two very different operations. 
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Modeled Turbulent Terms - LES 




PiMiUj 



Subgrid Stress 


p( Uj h - Ujh) - qj 5 


Subgrid heat flux 


These terms replace only turbulent effects that are smaller than the 
numerical scheme and grid (hence called subgrid) can resolve. 
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Comparison 


of RANS and LES Equations 
1 - Continuity w 



recHHO- 



+ ^-{pu,) = 0 
ax, j 

+ X"(^“') = 0 

ox i 


RANS 

LES 
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Comparison of RANS and LES Equations 

2 - Momentum 
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RANS 

LES 
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Comparison of RANS and LES Equations 

3 - Energy 



cE. d 
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* - A -T\ 3 /— . -I 

a acr J '-' , v ac^ J ' 1] 
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dE. d 
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3c 


LES 
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Major Differences Between Running 

RANS and LES 






All turbulent stresses are replaced by 
the averaged effect (numerically 
diffusive) 

Dominant turbulent stresses are 
calculated (numerically resolve 
unsteady behavior) - can't have too 
much diffusion or turbulence goes 
away. 

Constant CFL - goal is to get to 
convergence as fast as possible. 

Must run time accurately to capture 
time-varying turbulence. 

Grids are packed to regions of high 
mean shear (stretching OK) 

Best grids are uniform, isotropic - 
need to be of size to capture scales 
of interest. 

Numerical scheme designed for 
reasonable accuracy, shock wave 
capturing, convergence 
characteristics. 

Numerical scheme driven by need 
for high order of accuracy for 
resolving unsteady behavior. 
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RANS Turbulence Modeling 



• Reynolds-Averaged Navier-Stokes (RANS) - replaces all 
unsteady turbulent motion with modeled turbulent stresses. 

• Practical State of the art is two-equation models: k-e , k-co ,k-( 
Menter Shear-Stress Transport (SST) is popular “hybrid model” 
combining k-e and k-co. 

• For subsonic/transonic external aerodynamics, one equation 
models such as Spalart-Allmaras are popular - not used as 
much in propulsion flows. 

• Full Reynolds-Stress Models - offer more complete 
representation of 3-D turbulent stress field, but have not lived up 
to promise in terms of improved predictions. 

• Explicit algebraic stress models (EASMs) solve 2-eqn models, 
but used additional relations to obtain “Reynolds-stress-like” 
behavior. 
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Commonly Used RANS Models 






TECHKC*' 


1. Zero-equation (Algebraic): Cebeci-Smith, Baldwin-Lomax 

2. “Half-equation”: Johnson-King 

3. One-equation: Baldwin-Barth, Spalart-Allmaras 

4. Two-equation: 

a) k-8 : Jones-Launder (standard), Chien, many others 

b) k-co: Wilcox, Menter (BSL and SST) 

5. Explicit algebraic stress models (EASM) 

6. Reynolds-Stress closures 


The first 4 categories are “eddy viscosity models” where: 

f - - \ 


zl = -pu\u) = li 


2 ^_ 2 % 




J 


3 3c 


j 


J 
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Spalart-Allmaras One Equation RANS Model 
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Jones-Launder Two-Equation k-s RANS Model 
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• Jones-Launder form is referred to as the “standard” k-e model. 

• Works well for attached boundary layers. Underestimates size of flow 
separation (early reattachments). 

• Works reasonably well for mixing layers, jet flows. Compressibility 
corrections sometimes used for mixing problems at high convective 
Mach numbers. 



Menter Two-Equation k-co “Shear Stress 
Transport” (SST) RANS Model (1 of 2) 
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Menter Two-Equation k-co “Shear Stress 
Transport” (SST) RANS Model (2 of 2) 




TECH*®-" 


• SST term comes from eddy viscosity expression. The hybrid k-co, 
k-e model without this term is referred to as the “Baseline” or 
BSL model. 

• Inner model reduces to original Wilcox k-co formulation. Outer 
model comes from transformed k-e model and is supposed to 
reduce to the “standard” k-e model but differs in cross diffusion 
term and diffusion coefficients. 

• Model works quite well for attached boundary layers, mild 
separations, mixing flows (including jets). Numerically stable. 

• All k-co models may have an issue with artificial turbulence 
production behind strong normal shock. 
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Direct Calculation Methods 
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• Direct Numerical Simulation (DNS) - calculate all turbulent 
scales down to the Kolmogorov scale - impractical for 
engineering flows. 

• Large-Eddy Simulation (LES) - directly calculate largest scales 
and reserve modeling for smallest “subgrid-scale” stresses - 
active research showing promise in combustor and jet plume 
regions. 

• Hybrid RANS/LES - has become popular in recent years - most 
effective use has been for flows where RANS can be used in 
attached boundary layers and LES away from walls. 

- Demarcated or zonal hybrid RANS/LES - clear distinction is made between RANS 
and LES regions. Some physical mechanism is responsible for transition to 
turbulence. This was intent behind design of Detached Eddy Simulation (DES). 

- Continuous modeling - RANS and LES regions are not clearly separated - 
solution is expected to adjust, based on resolution. Desirable in theory, but difficult 
to achieve due to competing natures of RANS and LES. 
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Smagorinsky SGS Model 
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« M »« >4 
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f ^*Tii5S&- 0 


rf = ucAfp^Cs, vy - |c,a 2 p ^ 5 

p ss ‘ = p(CA) 2 pj7 J ' ! 

* = 

A = (AjcAyAz ) 1 ' 3 

• A few possibilities A = max(Ax>A) , jAz) 
for the subgrid 

(Ax)" + (Ay ) 2 + (Az ) 2 
3 

• Note similarity of functional form to mixing-length RANS model 
(i.e. Cebeci-Smith); gradient-diffusion formulation; eddy 
viscosity that adds to laminar viscosity just as is done in RANS. 

• The effect on N-S equations, however, is very different - only 
replacing subgrid turbulent stresses. 


turbulent length 

scale: A = 
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Transition Modeling 
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• Several RANS-based models tried over the past several years - some solving 
additional transport equations for intermittency, Re e . 

• Some success for flows with high freestream turbulence intensity - i.e. turbine 
cascades where bypass transition is dominant mechanism. 

• Modal growth situations not easily represented by RANS-based techniques. 

• Work shown here is with a model based on the Menter SST k-co turbulence model, 
with transition modifications by Langtry, Sjolander, & Menter. 

• Our work with the baseline published model indicated difficulties: (1) inability to 
reproduce experimentally observed transition, (2) significant grid sensitivity, (3) 
inability to become fully turbulent beyond transition. New formulation described in 
Denissen, Yoder, Georgiadis, NASA TM 2008-215451. 


TKE equation 
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dt 
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Modified model 
formulation: 
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J ;(3.28£ - 4 )Re v - (3.94£ - 7)Re‘i (1.43£ - 10)i?4] ; Re v < 1000 
| [0.12 4- {IE - 5 )Re v ] ; > 1000 
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Boundary Layer Transition Model 
Incompressible Validation 




Incompressible Validation: 

• Transition locations and skin friction examined for T3A 
benchmark data (ERCOFTAC) 

• Several freestream intensities investigated. 

• Grid sensitivity is high for incompressible cases. 


C f for FSTI = 2% 



C f Variation with FSTI 
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Boundary Layer Transition Model ^RKi 
Hypersonic Validation 


Hypersonic Validation: 

• Mach 7.93, 7 degree straight cone investigated in AEDC Tunnel B, T w / T 0 = 0.42. 

• Heat transfer measurements by Kimmel, JFE 1997. 

• Integrated heat transfer: Transition-SST (6.7% error), Fully turbulent SST (18.5 % 

error). i 1 
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Transition Modeling Conclusions 




• RANS-based models only applicable for bypass transition 
situations. 

• Free-flight transition is normally modal growth - a reliable RANS- 
based method is not likely promising. 

• LES is not promising either because accurately capturing the 
small disturbances is crucial - which LES will model/smear. 

• Long Term Prospects - DNS, e N methods. 
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Shock-Wave Turbulent Boundary 
Layer Interactions (SWTBLIs) 


Ihsst, 



• Pervasive to the entire hypersonic propulsion flowpath. 

• Major challenge to RANS, LES and hybrid RANS-LES techniques. 

• Nominally 2D problems are inherently 3D. 
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UFAST - Mach 2.25 Test Case 




/ HSST 

* .*-< *>« >■* t 

HKM S^CSO 
k SYSTEMS TEST * 


f l TECW* 0 - 


• 2010 AIAA Workshop: UFAST and U. of Michigan cases, targeted 
at representing supersonic aircraft inlets. 

• Several organizations submitted results - RANS, LES, hybrids 


Z 



I I 1 I 1 I i i I I 1 1 1 I i I i i i i I i i i i 1 i i i i 1 
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UFAST Velocity Contours 



Exploring minor change to Menter SST model’s stress limiter. 

Experiment Menter SST k-co, a 1 = 0.34 



Menter SST k-co 


Menter SST k-co, a 1 = 0.355 
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Menter SST k-co, a 1 = 0.37 
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Hypersonic International Flight Research (HIFiRE) ^ 
Direct-Connect Rig (HDCR) - Mach 5.8 Flight Case 


fHSST 

< M + < ¥•* I 

v K*M SPICO 




Mach No.: 0 0.5 1 1.5 2 2.5 


Expt, Body -Side Wall 
Expt, Cowl-Side Wall 
SST 

SST.al =0.355 














SWTBLI Modeling Conclusions 



e k-e models are generally overly optimistic on boundary layer 
health - smaller separations than expt. 

• k-t£> models usually work better for mild adverse pressure 
gradients, small separations, Menter SST predicts larger 
separations than expt. 

• One equation models (i.e. SA) provide similar accuracy to multi- 
equation models. 

• EASMs offer minimal improvement. 

• Some success using LES at AIAA Workshop, inflow conditions & 
matching Re are significant challenges. 

• Hybrid RANS-LES also being investigated - however, where is 
the switch from RANS to LES done? 
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Combustor/Exhaust System Modeling 
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SYSTEMS TEST 




• Several interacting phenomena - kinetics, turbulence, heat 
transfer, thermal-structural effects. 

• Practical state-of-the-art: Arrhenius form for reaction rates, 2 eqn 
turbulence model, constant Pr t , Sc t . Specified wall temperatures 
or heat fluxes. 

• Most practical scramjet experiments: only centerline pressures 
available; More data and/or unit problems are desirable. 

University of Virginia Supersonic 

Combustion Facility (UVA SCF): 

• Mach 5 enthalpy, Mach 2 isolator 

• overall pressure ratio - 4 

• H 2 fueled, clean air and vitiated air. 

• Documented heat transfer rates and 

wall temperatures. 

• NASA-sponsored experiments 

focused on mode transition behavior. 

• Continuing experiments through 

National Center. 

x/H : 
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Turbulent transport in energy 
and species equations 

nrn 

Turbulent heat flux: q T = — pu\h = —k 1 

<3c ; - 


Turbulent Prandtl number: 



Turbulent species flux: 


m] = -pu'.w j = -D n 


T 




T H 

Turbulent Schmidt number: Sc — — — 

d t 

The turbulent Prandtl and Schmidt number are frequently set 
equal to 0.9. However, it is believed that realistic values can 
be significantly different for many flows - particularly in 
extreme environments such as scramjets. 


40 

















Pr t and Sc t Sensitivity for USAF Scramjet 



/ USST t 




Fuel Injectors 

Cavity Flameholder 






Facility Nozzle 

Isolator 

Combustor 

Aft Nozzle 




An "optimized "Pr t and Sc t for one case do 

not guarantee optimal performance for Figure: courtesy of Robert A. Baurle, NASA LaRC 

other (p’s, turb. models, kinetics, etc. 
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Pr t Sensitivity for USAF Scramjet 
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>•,= 1.8 


Figure: Courtesy of Robert A. Baurle, NASA LaRC 




Burrows-Kurkov “Unit” Test Case 
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• Mach 2.4 vitiated air / sonic hydrogen experiment (1973). 

• Used extensively for investigations/validation of H 2 -air CFD methods 
(kinetics, variable Pr t , Sc t , hybrid RANS-LES...), perhaps overused. 

• Measurements of species concentrations and temperatures. 


Test section initial 


Test section intermediate 
measurement station, 
18. 3 cm 


Test section exit 
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Sc t Effects on Ignition Point for 
Burrow-Kurkov Test Case 






4 TECH* 0 ’' 0 


Pr t = 0.7 (constant) for all cases 


Sc t - 0.5 §. 

>» 


Sc t = 0.7 
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Hybrid RANS/LES Calculations of UVA 
Dual-Mode Scramjet, O = 0.17 



Temperature 



Eddy 

viscosity 



Figure: Courtesy of Jack Edwards, NCSU 
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Hybrid RANS/LES Calculations of UVA 
Dual-Mode Sc ramjet, d> = 0.17 


CARS comparisons (temperature): (X/H = 6, 12, 18) 


X/H=12 


X/H=18 


RANS 


LES/RANS 


LES/RANS 


Figure: Courtesy of Jack Edwards, NCSU 


(interpolated) 






Compressible Mixing 




• Most recent free shear layer mixing research has been in support of jet 
aeroacoustics research (subsonics and supersonics). 

• Practical state-of-the-art for RANS is also two-equation modeling. 

• Some research in variable Pr t for hot jet cases. 

• Most research support is towards LES-based methods. 

• Key LES issues: 

1 . Inflow boundary treatment 

2. Grid resolution/sensitivity 

3. Farfield noise propagation techniques. 
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Jets and Mixing - RANS 


0 - 4 ) 


Hssr 



-1 \ D Tt C\ v 


"«■ ( rtcwtO' 


RANS Findings: 

• RANS underpredict mixing for incompressible jets - initial shear layer is difficulty. 

• Uncorrected RANS models overpredict mixing rate for supersonic jets and mixing layers. 

• Effects of temperature and 3D jet effects are not modeled correctly. 

• Compressibility corrections (i.e. Sarkar) are highly empirical and do not reproduce correct 
fluid dynamic effects. 


Mach 0.5 Jet 


Mach 2.0 Jet 
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Jet Mixing - LES 


.<0 I /. 








in Him 
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• Acoustic Reference Nozzle (ARN) and Simple Metal Chevron (SMC) 
configurations - tested at GRC, investigated by several LES researchers. 

• Two Mach 0.9 jet simulations considered here: (1 ) DeBonis (GRC) DRP with 4 
stage RK, 3.5 - 9.2 million points and (2) Uzun (FSU), 4 th order compact scheme 
with 4 stage RK, 50 - 400 million points. 


DeBonis (GRC) grid: 










X/D : 
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Mean Axial Velocity 



srline Statistics (GRC) 






Axial Turbulent Intensity Radial Turbulent Intensity 



x/Dj x/Dj 


Figure: Courtesy of Jim DeBonis, NASA GRC 
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Turbulence Intensity Comparisons 
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SMC - Jet Decay and Acoustic Radiation 




Figure: Courtesy of AH Uzun, FSU 
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SMC - Jet Decay and Acoustic Radiation 


Figure: Courtesy of Ali Uzun, FSU 




Combustor/Exhaust System 
Modeling Enhancement Needs 





• RANS: 

- Better prediction of 3D, compressible mixing; highly separated/recirculating 
flow in flameholder/cavity, SWTBLIs, turbulent-chemistry interactions. 

- More accurate boundary conditions for thermal state. 

- Variable Pr t and Sc, capability. 

• LES: 

- Capability to handle wall bounded and free shear layer regions. Hybrid 
RANS/LES methods are under investigation - but location of RANS-to-LES 
switch has significant effect. 

- Significant uncertainty remains in how to best perform jet/mixing 
simulations. Highly desirable to establish “best practices” if possible. 

- Models for turbulent/chemistry interactions, i.e. Filtered Density Functions 
(FDFs). 
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Experimental - Validation Data Needs 




• Centerline pressure distributions are not sufficient for validation / 
calibration of turbulent flow CFD. There are too many interacting 
features in scramjet flowpaths - unlike subsonic/transonic 
aerodynamics. 

• More complete turbulent statistics for momentum, thermal, and 
species transport are needed. 

• Advanced Diagnostics: CARS, PLIF, PIV - for unit problems, then 
more complex cases. 

• Quantify uncertainty - e.g. PIV is powerful technique, but prone to 
high uncertainty in crucial regions such as initial mixing regions. 

• Consider revisiting experiments such as Burrows-Kurkov with the 
advanced techniques. 

• Design experiments to avoid contamination of focus region - i.e. 
SWBLI cases - nearly all experiments are in small tunnels where 
sidewall separations dominate region of interest. 
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Conclusions 



• Many extremely difficult challenges remain in turbulence modeling 
for air-breathing propulsion flows. 


• Status of RANS Modeling for high speed propulsion flowpaths: 
Not much advancement in practical state-of-the-art in 2 decades. 

• Dominant features of 3-D flow, large separations, SWTBLIs, 
chemically reacting flow, compressibility, turbulent transport of 
heat and species - overwhelm the capabilities of current RANS 
methods. 

• Tweaking one turbulence modeling parameter while holding all 
others fixed until centerline pressure distribution matches 

experimental data (typical practice for scramjets) is of minimal 
value. 

• LES and related methods are demonstrating some promise, but 
have their own modeling issues and (1) are not of sufficient 
maturity for most problems, (2) computing power is not readily 
available to use in a production engineering environment, (3) 
minimal consistency between groups in how to achieve most 
accurate results. 
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